A framework for conducting time-varying genome-wide association studies: An application to body mass index across childhood in six multiethnic cohorts.

Genetic effects on changes in human traits over time are understudied and may have important pathophysiological impact. We propose a framework that enables data quality control, implements mixed models to evaluate trajectories of change in traits, and estimates phenotypes to identify age-varying genetic effects in genome-wide association studies (GWASs). Using childhood body mass index (BMI) as an example, we included 71,336 participants from six cohorts and estimated the slope and area under the BMI curve within four time periods (infancy, early childhood, late childhood and adolescence) for each participant, in addition to the age and BMI at the adiposity peak and the adiposity rebound. GWAS on each of the estimated phenotypes identified 28 genome-wide significant variants at 13 loci across the 12 estimated phenotypes, one of which was novel (in DAOA) and had not been previously associated with childhood or adult BMI. Genetic studies of changes in human traits over time could uncover novel biological mechanisms influencing quantitative traits.

Pathogenic, Total Loss-of-Function DYRK1B Variants Cause Monogenic Obesity Associated With Type 2 Diabetes.

OBJECTIVE: Rare variants in DYRK1B have been described in some patients with central obesity, type 2 diabetes, and early-onset coronary disease. Owing to the limited number of conducted studies, the broader impact of DYRK1B variants on a larger scale has yet to be investigated. RESEARCH DESIGN AND METHODS: DYRK1B was sequenced in 9,353 participants from a case-control study for obesity and type 2 diabetes. Each DYRK1B variant was functionally assessed in vitro. Variant pathogenicity was determined using criteria from the American College of Medical Genetics and Genomics (ACMG). The effect of pathogenic or likely pathogenic (P/LP) variants on metabolic traits was assessed using adjusted mixed-effects score tests. RESULTS: Sixty-five rare, heterozygous DYRK1B variants were identified and were not associated with obesity or type 2 diabetes. Following functional analyses, 20 P/LP variants were pinpointed, including 6 variants that exhibited a fully inhibitory effect (P/LP-null) on DYRK1B activity. P/LP and P/LP-null DYRK1B variants were associated with increased BMI and obesity risk; however, the impact was notably more pronounced for the P/LP-null variants (effect of 8.0 ± 3.2 and odds ratio of 7.9 [95% CI 1.2-155]). Furthermore, P/LP-null variants were associated with higher fasting glucose and type 2 diabetes risk (effect of 2.9 ± 1.0 and odds ratio of 4.8 [95% CI 0.85-37]), while P/LP variants had no effect on glucose homeostasis. CONCLUSIONS: P/LP, total loss-of-function DYRK1B variants cause monogenic obesity associated with type 2 diabetes. This study underscores the significance of conducting functional assessments in order to accurately ascertain the tangible effects of P/LP DYRK1B variants.

Pathogenic monoallelic variants in GLIS3 increase type 2 diabetes risk and identify a subgroup of patients sensitive to sulfonylureas.

AIMS/HYPOTHESIS: GLIS3 encodes a transcription factor involved in pancreatic beta cell development and function. Rare pathogenic, bi-allelic mutations in GLIS3 cause syndromic neonatal diabetes whereas frequent SNPs at this locus associate with common type 2 diabetes risk. Because rare, functional variants located in other susceptibility genes for type 2 diabetes have already been shown to strongly increase individual risk for common type 2 diabetes, we aimed to investigate the contribution of rare pathogenic GLIS3 variants to type 2 diabetes. METHODS: GLIS3 was sequenced in 5471 individuals from the Rare Variants Involved in Diabetes and Obesity (RaDiO) study. Variant pathogenicity was assessed following the criteria established by the American College of Medical Genetics and Genomics (ACMG). To address the pathogenic strong criterion number 3 (PS3), we conducted functional investigations of these variants using luciferase assays, focusing on capacity of GLIS family zinc finger 3 (GLIS3) to bind to and activate the INS promoter. The association between rare pathogenic or likely pathogenic (P/LP) variants and type 2 diabetes risk (and other metabolic traits) was then evaluated. A meta-analysis combining association results from RaDiO, the 52K study (43,125 individuals) and the TOPMed study (44,083 individuals) was finally performed. RESULTS: Through targeted resequencing of GLIS3, we identified 105 rare variants that were carried by 395 participants from RaDiO. Among them, 49 variants decreased the activation of the INS promoter. Following ACMG criteria, 18 rare variants were classified as P/LP, showing an enrichment in the last two exons compared with the remaining exons (p<5×10-6; OR>3.5). The burden of these P/LP variants was strongly higher in individuals with type 2 diabetes (p=3.0×10-3; OR 3.9 [95% CI 1.4, 12]), whereas adiposity, age at type 2 diabetes diagnosis and cholesterol levels were similar between variant carriers and non-carriers with type 2 diabetes. Interestingly, all carriers with type 2 diabetes were sensitive to oral sulfonylureas. A total of 7 P/LP variants were identified in both 52K and TOPMed studies. The meta-analysis of association studies obtained from RaDiO, 52K and TOPMed showed an enrichment of P/LP GLIS3 variants in individuals with type 2 diabetes (p=5.6×10-5; OR 2.1 [95% CI 1.4, 2.9]). CONCLUSIONS/INTERPRETATION: Rare P/LP GLIS3 variants do contribute to type 2 diabetes risk. The variants located in the distal part of the protein could have a direct effect on its functional activity by impacting its transactivation domain, by homology with the mouse GLIS3 protein. Furthermore, rare P/LP GLIS3 variants seem to have a direct clinical effect on beta cell function, which could be improved by increasing insulin secretion via the use of sulfonylureas.

Exploring the role of purinergic receptor P2RY1 in type 2 diabetes risk and pathophysiology: Insights from human functional genomics.

OBJECTIVE: Human functional genomics has proven powerful in discovering drug targets for common metabolic disorders. Through this approach, we investigated the involvement of the purinergic receptor P2RY1 in type 2 diabetes (T2D). METHODS: P2RY1 was sequenced in 9,266 participants including 4,177 patients with T2D. In vitro analyses were then performed to assess the functional effect of each variant. Expression quantitative trait loci (eQTL) analysis was performed in pancreatic islets from 103 pancreatectomized individuals. The effect of P2RY1 on glucose-stimulated insulin secretion was finally assessed in human pancreatic beta cells (EndoCßH5), and RNA sequencing was performed on these cells. RESULTS: Sequencing P2YR1 in 9,266 participants revealed 22 rare variants, seven of which were loss-of-function according to our in vitro analyses. Carriers, except one, exhibited impaired glucose control. Our eQTL analysis of human islets identified P2RY1 variants, in a beta-cell enhancer, linked to increased P2RY1 expression and reduced T2D risk, contrasting with variants located in a silent region associated with decreased P2RY1 expression and increased T2D risk. Additionally, a P2RY1-specific agonist increased insulin secretion upon glucose stimulation, while the antagonist led to decreased insulin secretion. RNA-seq highlighted TXNIP as one of the main transcriptomic markers of insulin secretion triggered by P2RY1 agonist. CONCLUSION: Our findings suggest that P2RY1 inherited or acquired dysfunction increases T2D risk and that P2RY1 activation stimulates insulin secretion. Selective P2RY1 agonists, impermeable to the blood-brain barrier, could serve as potential insulin secretagogues.

Gene-educational attainment interactions in a multi-population genome-wide meta-analysis identify novel lipid loci.

Introduction: Educational attainment, widely used in epidemiologic studies as a surrogate for socioeconomic status, is a predictor of cardiovascular health outcomes. Methods: A two-stage genome-wide meta-analysis of low-density lipoprotein cholesterol (LDL), high-density lipoprotein cholesterol (HDL), and triglyceride (TG) levels was performed while accounting for gene-educational attainment interactions in up to 226,315 individuals from five population groups. We considered two educational attainment variables: "Some College" (yes/no, for any education beyond high school) and "Graduated College" (yes/no, for completing a 4-year college degree). Genome-wide significant (p < 5 × 10-8) and suggestive (p < 1 × 10-6) variants were identified in Stage 1 (in up to 108,784 individuals) through genome-wide analysis, and those variants were followed up in Stage 2 studies (in up to 117,531 individuals). Results: In combined analysis of Stages 1 and 2, we identified 18 novel lipid loci (nine for LDL, seven for HDL, and two for TG) by two degree-of-freedom (2 DF) joint tests of main and interaction effects. Four loci showed significant interaction with educational attainment. Two loci were significant only in cross-population analyses. Several loci include genes with known or suggested roles in adipose (FOXP1, MBOAT4, SKP2, STIM1, STX4), brain (BRI3, FILIP1, FOXP1, LINC00290, LMTK2, MBOAT4, MYO6, SENP6, SRGAP3, STIM1, TMEM167A, TMEM30A), and liver (BRI3, FOXP1) biology, highlighting the potential importance of brain-adipose-liver communication in the regulation of lipid metabolism. An investigation of the potential druggability of genes in identified loci resulted in five gene targets shown to interact with drugs approved by the Food and Drug Administration, including genes with roles in adipose and brain tissue. Discussion: Genome-wide interaction analysis of educational attainment identified novel lipid loci not previously detected by analyses limited to main genetic effects.

Human GLP1R variants affecting GLP1R cell surface expression are associated with impaired glucose control and increased adiposity.

The glucagon-like peptide 1 receptor (GLP1R) is a major drug target with several agonists being prescribed in individuals with type 2 diabetes and obesity1,2. The impact of genetic variability of GLP1R on receptor function and its association with metabolic traits are unclear with conflicting reports. Here, we show an unexpected diversity of phenotypes ranging from defective cell surface expression to complete or pathway-specific gain of function (GoF) and loss of function (LoF), after performing a functional profiling of 60 GLP1R variants across four signalling pathways. The defective insulin secretion of GLP1R LoF variants is rescued by allosteric GLP1R ligands or high concentrations of exendin-4/semaglutide in INS-1 823/3 cells. Genetic association studies in 200,000 participants from the UK Biobank show that impaired GLP1R cell surface expression contributes to poor glucose control and increased adiposity with increased glycated haemoglobin A1c and body mass index. This study defines impaired GLP1R cell surface expression as a risk factor for traits associated with type 2 diabetes and obesity and provides potential treatment options for GLP1R LoF variant carriers.

Bidirectional Mendelian Randomization and Multiphenotype GWAS Show Causality and Shared Pathophysiology Between Depression and Type 2 Diabetes.

OBJECTIVE: Depression is a common comorbidity of type 2 diabetes. We assessed the causal relationships and shared genetics between them. RESEARCH DESIGN AND METHODS: We applied two-sample, bidirectional Mendelian randomization (MR) to assess causality between type 2 diabetes and depression. We investigated potential mediation using two-step MR. To identify shared genetics, we performed 1) genome-wide association studies (GWAS) separately and 2) multiphenotype GWAS (MP-GWAS) of type 2 diabetes (19,344 case subjects, 463,641 control subjects) and depression using major depressive disorder (MDD) (5,262 case subjects, 86,275 control subjects) and self-reported depressive symptoms (n = 153,079) in the UK Biobank. We analyzed expression quantitative trait locus (eQTL) data from public databases to identify target genes in relevant tissues. RESULTS: MR demonstrated a significant causal effect of depression on type 2 diabetes (odds ratio 1.26 [95% CI 1.11-1.44], P = 5.46 × 10-4) but not in the reverse direction. Mediation analysis indicated that 36.5% (12.4-57.6%, P = 0.0499) of the effect from depression on type 2 diabetes was mediated by BMI. GWAS of type 2 diabetes and depressive symptoms did not identify shared loci. MP-GWAS identified seven shared loci mapped to TCF7L2, CDKAL1, IGF2BP2, SPRY2, CCND2-AS1, IRS1, CDKN2B-AS1. MDD has not brought any significant association in either GWAS or MP-GWAS. Most MP-GWAS loci had an eQTL, including single nucleotide polymorphisms implicating the cell cycle gene CCND2 in pancreatic islets and brain and the insulin signaling gene IRS1 in adipose tissue, suggesting a multitissue and pleiotropic underlying mechanism. CONCLUSIONS: Our results highlight the importance to prevent type 2 diabetes at the onset of depressive symptoms and the need to maintain a healthy weight in the context of its effect on depression and type 2 diabetes comorbidity.

Heterozygous pathogenic variants in POMC are not responsible for monogenic obesity: Implication for MC4R agonist use.

PURPOSE: Recessive deficiency of proopiomelanocortin (POMC) causes childhood-onset severe obesity. Cases can now benefit from the melanocortin 4 receptor agonist setmelanotide. Furthermore, a phase 3 clinical trial is evaluating setmelanotide in heterozygotes for POMC. We performed a large-scale genetic analysis to assess the effect of heterozygous, pathogenic POMC variants on obesity. METHODS: A genetic analysis was performed in a family including 2 cousins with childhood-onset obesity. We analyzed the obesity status of heterozygotes for pathogenic POMC variants in the Human Gene Mutation Database. The association between heterozygous pathogenic POMC variants and obesity risk was assessed using 190,000 exome samples from UK Biobank. RESULTS: The 2 cousins carried a compound heterozygous pathogenic variant in POMC. Six siblings were heterozygotes; only 1 of them had obesity. In Human Gene Mutation Database, we identified 60 heterozygotes for pathogenic POMC variants, of whom 14 had obesity. In UK Biobank, heterozygous pathogenic POMC variants were not associated with obesity risk, but they modestly increased body mass index levels. CONCLUSION: Heterozygous pathogenic POMC variants do not contribute to monogenic obesity, but they slightly increase body mass index. Setmelanotide use in patients with obesity, which would only be based on the presence of a heterozygous POMC variant, can be questioned.

Biallelic Mutations in P4HTM Cause Syndromic Obesity.

We previously demonstrated that 50% of children with obesity from consanguineous families from Pakistan carry pathogenic variants in known monogenic obesity genes. Here, we have discovered a novel monogenetic recessive form of severe childhood obesity using an in-house computational staged approach. The analysis included whole-exome sequencing data of 366 children with severe obesity, 1,000 individuals of the Pakistan Risk of Myocardial Infarction Study (PROMIS) study, and 200,000 participants of the UK Biobank to prioritize genes harboring rare homozygous variants with putative effect on human obesity. We identified five rare or novel homozygous missense mutations predicted deleterious in five consanguineous families in P4HTM encoding prolyl 4-hydroxylase transmembrane (P4H-TM). We further found two additional homozygous missense mutations in children with severe obesity of Indian and Moroccan origin. Molecular dynamics simulation suggested that these mutations destabilized the active conformation of the substrate binding domain. Most carriers also presented with hypotonia, cognitive impairment, and/or developmental delay. Three of the five probands died of pneumonia during the first 2 years of the follow-up. P4HTM deficiency is a novel form of syndromic obesity, affecting 1.5% of our children with obesity associated with high mortality. P4H-TM is a hypoxia-inducible factor that is necessary for survival and adaptation under oxygen deprivation, but the role of this pathway in energy homeostasis and obesity pathophysiology remains to be elucidated.

Contribution of heterozygous PCSK1 variants to obesity and implications for precision medicine: a case-control study.

BACKGROUND: Rare biallelic pathogenic mutations in PCSK1 (encoding proprotein convertase subtilisin/kexin type 1 [PC1/3]) cause early-onset obesity associated with various endocrinopathies. Setmelanotide has been approved for carriers of these biallelic mutations in the past 3 years. We aimed to perform a large-scale functional genomic study focusing on rare heterozygous variants of PCSK1 to decipher their putative impact on obesity risk. METHODS: This case-control study included all participants with overweight and obesity (ie, cases) or healthy weight (ie, controls) from the RaDiO study of three community-based and one hospital-based cohort in France recruited between Jan 1, 1995, and Dec 31, 2000. In adults older than 18 years, healthy weight was defined as BMI of less than 25·0 kg/m2, overweight as 25·0-29·9 kg/m2, and obesity as 30·0 kg/m2 or higher. Participants with type 2 diabetes had fasting glucose of 7·0 mmol/L or higher or used treatment for hyperglycaemia (or both) and were negative for islet or insulin autoantibodies. Functional assessment of rare missense variants of PCSK1 was performed. Pathogenicity clusters of variants were determined with machine learning. The effect of each cluster of PCSK1 variants on obesity was assessed using the adjusted mixed-effects score test. FINDINGS: All 13 coding exons of PCSK1 were sequenced in 9320 participants (including 7260 adults and 2060 children and adolescents) recruited from the RaDiO study. We detected 65 rare heterozygous PCSK1 variants, including four null variants and 61 missense variants that were analysed in vitro and clustered into five groups (A-E), according to enzymatic activity. Compared with the wild-type, 15 missense variants led to complete PC1/3 loss of function (group A; reference) and rare exome variant ensemble learner (REVEL) led to 15 (25%) false positives and four (7%) false negatives. Carrying complete loss-of-function or null PCSK1 variants was significantly associated with obesity (six [86%] of seven carriers vs 1518 [35%] of 4395 non-carriers; OR 9·3 [95% CI 1·5-177·4]; p=0·014) and higher BMI (32·0 kg/m2 [SD 9·3] in carriers vs 27·3 kg/m2 [6·5] in non-carriers; mean effect π 6·94 [SE 1·95]; p=0·00029). Clusters of PCSK1 variants with partial or neutral effect on PC1/3 activity did not have an effect on obesity or overweight and on BMI. INTERPRETATION: Only carriers of heterozygous, null, or complete loss-of-function PCSK1 variants cause monogenic obesity and, therefore, might be eligible for setmelanotide. In silico tests were unable to accurately detect these variants, which suggests that in vitro assays are necessary to determine the variant pathogenicity for genetic diagnosis and precision medicine purposes. FUNDING: Agence Nationale de la Recherche, European Research Council, National Center for Precision Diabetic Medicine, European Regional Development Fund, Hauts-de-France Regional Council, and the European Metropolis of Lille.

Rare Variant Analysis of Obesity-Associated Genes in Young Adults With Severe Obesity From a Consanguineous Population of Pakistan.

Recent advances in genetic analysis have significantly helped in progressively attenuating the heritability gap of obesity and have brought into focus monogenic variants that disrupt the melanocortin signaling. In a previous study, next-generation sequencing revealed a monogenic etiology in ∼50% of the children with severe obesity from a consanguineous population in Pakistan. Here we assess rare variants in obesity-causing genes in young adults with severe obesity from the same region. Genomic DNA from 126 randomly selected young adult obese subjects (BMI 37.2 ± 0.3 kg/m2; age 18.4 ± 0.3 years) was screened by conventional or augmented whole-exome analysis for point mutations and copy number variants (CNVs). Leptin, insulin, and cortisol levels were measured by ELISA. We identified 13 subjects carrying 13 different pathogenic or likely pathogenic variants in LEPR, PCSK1, MC4R, NTRK2, POMC, SH2B1, and SIM1. We also identified for the first time in the human, two homozygous stop-gain mutations in ASNSD1 and IFI16 genes. Inactivation of these genes in mouse models has been shown to result in obesity. Additionally, we describe nine homozygous mutations (seven missense, one stop-gain, and one stop-loss) and four copy-loss CNVs in genes or genomic regions previously linked to obesity-associated traits by genome-wide association studies. Unexpectedly, in contrast to obese children, pathogenic mutations in LEP and LEPR were either absent or rare in this cohort of young adults. High morbidity and mortality risks and social disadvantage of children with LEP or LEPR deficiency may in part explain this difference between the two cohorts.

Identification of Key Regions Mediating Human Melatonin Type 1 Receptor Functional Selectivity Revealed by Natural Variants.

Melatonin is a hormone mainly produced by the pineal gland and MT1 is one of the two G protein-coupled receptors (GPCRs) mediating its action. Despite an increasing number of available GPCR crystal structures, the molecular mechanism of activation of a large number of receptors, including MT1, remains poorly understood. The purpose of this study is to elucidate the structural elements involved in the process of MT1's activation using naturally occurring variants affecting its function. Thirty-six nonsynonymous variants, including 34 rare ones, were identified in MTNR1A (encoding MT1) from a cohort of 8687 individuals and their signaling profiles were characterized using Bioluminescence Resonance Energy Transfer-based sensors probing 11 different signaling pathways. Computational analysis of the experimental data allowed us to group the variants in clusters according to their signaling profiles and to analyze the position of each variant in the context of the three-dimensional structure of MT1 to link functional selectivity to structure. MT1 variant signaling profiles revealed three clusters characterized by (1) wild-type-like variants, (2) variants with selective defect of ßarrestin-2 recruitment, and (3) severely defective variants on all pathways. Our structural analysis allows us to identify important regions for ßarrestin-2 recruitment as well as for Gα12 and Gα15 activation. In addition to identifying MT1 domains differentially controlling the activation of the various signaling effectors, this study illustrates how natural variants can be used as tools to study the molecular mechanisms of receptor activation.

Gene-educational attainment interactions in a multi-ancestry genome-wide meta-analysis identify novel blood pressure loci.

Educational attainment is widely used as a surrogate for socioeconomic status (SES). Low SES is a risk factor for hypertension and high blood pressure (BP). To identify novel BP loci, we performed multi-ancestry meta-analyses accounting for gene-educational attainment interactions using two variables, "Some College" (yes/no) and "Graduated College" (yes/no). Interactions were evaluated using both a 1 degree of freedom (DF) interaction term and a 2DF joint test of genetic and interaction effects. Analyses were performed for systolic BP, diastolic BP, mean arterial pressure, and pulse pressure. We pursued genome-wide interrogation in Stage 1 studies (N = 117 438) and follow-up on promising variants in Stage 2 studies (N = 293 787) in five ancestry groups. Through combined meta-analyses of Stages 1 and 2, we identified 84 known and 18 novel BP loci at genome-wide significance level (P < 5 × 10-8). Two novel loci were identified based on the 1DF test of interaction with educational attainment, while the remaining 16 loci were identified through the 2DF joint test of genetic and interaction effects. Ten novel loci were identified in individuals of African ancestry. Several novel loci show strong biological plausibility since they involve physiologic systems implicated in BP regulation. They include genes involved in the central nervous system-adrenal signaling axis (ZDHHC17, CADPS, PIK3C2G), vascular structure and function (GNB3, CDON), and renal function (HAS2 and HAS2-AS1, SLIT3). Collectively, these findings suggest a role of educational attainment or SES in further dissection of the genetic architecture of BP.

Monogenic diabetes characteristics in a transnational multicenter study from Mediterranean countries.

BACKGROUND: Diagnosis of monogenic diabetes has important clinical implications for treatment and health expenditure. However, its prevalence remains to be specified in many countries, particularly from South Europe, North Africa and Middle-East, where non-autoimmune diabetes in young adults is increasing dramatically. AIMS: To identify cases of monogenic diabetes in young adults from Mediterranean countries and assess the specificities between countries. METHODS: We conducted a transnational multicenter study based on exome sequencing in 204 unrelated patients with diabetes (age-at-diagnosis: 26.1 ± 9.1 years). Rare coding variants in 35 targeted genes were evaluated for pathogenicity. Data were analyzed using one-way ANOVA, chi-squared test and factor analysis of mixed data. RESULTS: Forty pathogenic or likely pathogenic variants, 14 of which novel, were identified in 36 patients yielding a genetic diagnosis rate of 17.6%. The majority of cases were due to GCK, HNF1A, ABCC8 and HNF4A variants. We observed highly variable diagnosis rates according to countries, with association to genetic ancestry. Lower body mass index and HbA1c at study inclusion, and less frequent insulin treatment were hallmarks of pathogenic variant carriers. Treatment changes following genetic diagnosis have been made in several patients. CONCLUSIONS: Our data from patients in several Mediterranean countries highlight a broad clinical and genetic spectrum of diabetes, showing the relevance of wide genetic testing for personalized care of early-onset diabetes.

Pathogenic variants in actionable MODY genes are associated with type 2 diabetes.

Genome-wide association studies have identified 240 independent loci associated with type 2 diabetes (T2D) risk, but this knowledge has not advanced precision medicine. In contrast, the genetic diagnosis of monogenic forms of diabetes (including maturity-onset diabetes of the young (MODY)) are textbook cases of genomic medicine. Recent studies trying to bridge the gap between monogenic diabetes and T2D have been inconclusive. Here, we show a significant burden of pathogenic variants in genes linked with monogenic diabetes among people with common T2D, particularly in actionable MODY genes, thus implying that there should be a substantial change in care for carriers with T2D. We show that, among 74,629 individuals, this burden is probably driven by the pathogenic variants found in GCK, and to a lesser extent in HNF4A, KCNJ11, HNF1B and ABCC8. The carriers with T2D are leaner, which evidences a functional metabolic effect of these mutations. Pathogenic variants in actionable MODY genes are more frequent than was previously expected in common T2D. These results open avenues for future interventions assessing the clinical interest of these pathogenic mutations in precision medicine.

General regression model: A "model-free" association test for quantitative traits allowing to test for the underlying genetic model.

Most genome-wide association studies used genetic-model-based tests assuming an additive mode of inheritance, leading to underpowered association tests in case of departure from additivity. The general regression model (GRM) association test proposed by Fisher and Wilson in 1980 makes no assumption on the genetic model. Interestingly, it also allows formal testing of the underlying genetic model. We conducted a simulation study of quantitative traits to compare the power of the GRM test to the classical linear regression tests, the maximum of the three statistics (MAX), and the allele-based (allelic) tests. Simulations were performed on two samples sizes, using a large panel of genetic models, varying genetic models, minor allele frequencies, and the percentage of explained variance. In case of departure from additivity, the GRM was more powerful than the additive regression tests (power gain reaching 80%) and had similar power when the true model is additive. GRM was also as or more powerful than the MAX or allelic tests. The true simulated model was mostly retained by the GRM test. Application of GRM to HbA1c illustrates its gain in power. To conclude, GRM increases power to detect association for quantitative traits, allows determining the genetic model and is easily applicable.

Loss-of-function mutations in MRAP2 are pathogenic in hyperphagic obesity with hyperglycemia and hypertension.

The G-protein-coupled receptor accessory protein MRAP2 is implicated in energy control in rodents, notably via the melanocortin-4 receptor1. Although some MRAP2 mutations have been described in people with obesity1-3, their functional consequences on adiposity remain elusive. Using large-scale sequencing of MRAP2 in 9,418 people, we identified 23 rare heterozygous variants associated with increased obesity risk in both adults and children. Functional assessment of each variant shows that loss-of-function MRAP2 variants are pathogenic for monogenic hyperphagic obesity, hyperglycemia and hypertension. This contrasts with other monogenic forms of obesity characterized by excessive hunger, including melanocortin-4 receptor deficiency, that present with low blood pressure and normal glucose tolerance4. The pleiotropic metabolic effect of loss-of-function mutations in MRAP2 might be due to the failure of different MRAP2-regulated G-protein-coupled receptors in various tissues including pancreatic islets.

A catalog of genetic loci associated with kidney function from analyses of a million individuals.

Chronic kidney disease (CKD) is responsible for a public health burden with multi-systemic complications. Through trans-ancestry meta-analysis of genome-wide association studies of estimated glomerular filtration rate (eGFR) and independent replication (n = 1,046,070), we identified 264 associated loci (166 new). Of these, 147 were likely to be relevant for kidney function on the basis of associations with the alternative kidney function marker blood urea nitrogen (n = 416,178). Pathway and enrichment analyses, including mouse models with renal phenotypes, support the kidney as the main target organ. A genetic risk score for lower eGFR was associated with clinically diagnosed CKD in 452,264 independent individuals. Colocalization analyses of associations with eGFR among 783,978 European-ancestry individuals and gene expression across 46 human tissues, including tubulo-interstitial and glomerular kidney compartments, identified 17 genes differentially expressed in kidney. Fine-mapping highlighted missense driver variants in 11 genes and kidney-specific regulatory variants. These results provide a comprehensive priority list of molecular targets for translational research.

A multi-ancestry genome-wide study incorporating gene-smoking interactions identifies multiple new loci for pulse pressure and mean arterial pressure.

Elevated blood pressure (BP), a leading cause of global morbidity and mortality, is influenced by both genetic and lifestyle factors. Cigarette smoking is one such lifestyle factor. Across five ancestries, we performed a genome-wide gene-smoking interaction study of mean arterial pressure (MAP) and pulse pressure (PP) in 129 913 individuals in stage 1 and follow-up analysis in 480 178 additional individuals in stage 2. We report here 136 loci significantly associated with MAP and/or PP. Of these, 61 were previously published through main-effect analysis of BP traits, 37 were recently reported by us for systolic BP and/or diastolic BP through gene-smoking interaction analysis and 38 were newly identified (P < 5 × 10-8, false discovery rate < 0.05). We also identified nine new signals near known loci. Of the 136 loci, 8 showed significant interaction with smoking status. They include CSMD1 previously reported for insulin resistance and BP in the spontaneously hypertensive rats. Many of the 38 new loci show biologic plausibility for a role in BP regulation. SLC26A7 encodes a chloride/bicarbonate exchanger expressed in the renal outer medullary collecting duct. AVPR1A is widely expressed, including in vascular smooth muscle cells, kidney, myocardium and brain. FHAD1 is a long non-coding RNA overexpressed in heart failure. TMEM51 was associated with contractile function in cardiomyocytes. CASP9 plays a central role in cardiomyocyte apoptosis. Identified only in African ancestry were 30 novel loci. Our findings highlight the value of multi-ancestry investigations, particularly in studies of interaction with lifestyle factors, where genomic and lifestyle differences may contribute to novel findings.

Multi-ancestry genome-wide gene-smoking interaction study of 387,272 individuals identifies new loci associated with serum lipids.

The concentrations of high- and low-density-lipoprotein cholesterol and triglycerides are influenced by smoking, but it is unknown whether genetic associations with lipids may be modified by smoking. We conducted a multi-ancestry genome-wide gene-smoking interaction study in 133,805 individuals with follow-up in an additional 253,467 individuals. Combined meta-analyses identified 13 new loci associated with lipids, some of which were detected only because association differed by smoking status. Additionally, we demonstrate the importance of including diverse populations, particularly in studies of interactions with lifestyle factors, where genomic and lifestyle differences by ancestry may contribute to novel findings.